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A PANORAMIC VIEW OF THE MILKY WAY ANALOGUE NGC 891 
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ABSTRACT 

Recent panoramic observations of the dominant spiral galaxies of the Local Group have revolutionized 
our view of how these galaxies assemble their mass. However, it remains completely unclear whether the 
properties of the outer regions of the Local Group large spirals are typical. Here, we present the first 
panoramic view of a spiral galaxy beyond the Local Group, based on the largest, contiguous, ground- 
based imaging survey to date resolving the stellar halo of the nearest prime analogue of the Milky Way, 
NGC 891 (D « 10 Mpc). The low surface brightness outskirts of this galaxy are populated by multiple, 
coherent, and vast substructures over the ~ 90 kpc x 90 kpc extent of the survey. These include a giant 
stream, the first to be resolved into stars beyond the Local Group using ground-based facilities, that 
loops around the parent galaxy up to distances of ~ 50 kpc. The bulge and the disk of the galaxy are 
found to be surrounded by a previously undetected large, flat and thick cocoon-like stellar structure at 
vertical and radial distances of up to ~ 15 kpc and ~ 40 kpc respectively. 

Subject headings: galaxies: halos - galaxies: stellar content - galaxies: individual (NGC 891) 



1. INTRODUCTION 

In recent years, it has been increasingly recognized that 
many of the clues to the problem o f galaxy formation are 
preserved in galaxy outskirts (e.g. iJohnston et al.l [2008. 
and references therein). The current consensus is that 
large spirals begin as small fluctuations in the early Uni- 
verse and gro w by in situ star for mation and hierarchical 
merging fe.g. IWhite fc Reedll978l) . Subsequently, once a 
spiral is the dominant component in such mergers, it con- 
tinues accreting and disrupting sub-halos falling into its 
potential well. With the accumulation of accretion and 
disruption ev ents, massive spira ls build up stellar and dark 
matter halos (Abadi ct al. 200(|). The complete disruption 
of the sub-halos may take several orbits, distributing thus 
the tidally stripped stars over a broad range of distances 
from the main galaxy. As galaxies are predicted to assem- 
ble their outskirts from the disruption of a large number of 
satellites, these regions are expected to possess signifi cant 
density and chemical substructures (jFont et al.ll2006l) . 

Observationally however, the nature and the origin of 
those regions remains elusive. Much of what we know 
about their properties is base d on observations of the Local 
Grou p massive spirals (e.g. iFreeman fc Bland-Hawthorn! 
2002, and refer ences therein). Recent surveys find evi- 
dence that the Galaxy stellar halo is divisible into two 
components, with a moderately flat inner regions show- 
ing a modest prograde rotation, whereas the outer regions 
are less chemically evolved and exhibit a nearly spheri- 
cal d istribution with a retrograde rotation (jCarollo et al.l 
2008). Wide-field imaging data indicates th at the stel- 
lar halo of the Galaxy i s highly structured (llbata et al.l 
l2003t lYannv et alj|2003t iBelokurov et~aT1l2007D . suggest- 
ing that a larg e fraction of the Halo has been accreted 
from satellites ([Bell et al.l [20081 ). The outer regions of 
Andromeda have been recently observed to contain even 
more substructure and streams than observed around the 



Milky Way (jlbata et al.ll2007D . suggesting that its accre- 
tion history may have been mor e active than the sus- 
pected quiet on e of the Milky Way (|Mouhcine et al.ll2005at 
Ham mer et al.l 120071 ). It is completely unknown however 
if these properties are generic features of the outskirts of 
spirals, or are reflecting peculiar assembly histories. 

We have seen recently a dramatic progress in large-scale 
mapping of the Milky Way and Andromeda, with a num- 
ber of large surveys measuring photometric, kinematic, 
and chemical properties of individual stars over a wide 
galactic volume. Although those surveys will significantly 
advance our understanding of the assembly of these galax- 
ies, it cannot be assumed that a sample of two galaxies will 
provide the definitive solution to the nature and origin of 
the stellar content in the outskirts of spirals. The funda- 
mental next step to fully exploit the Local Group surveys 
and thereby to establish a comprehensive picture of the 
assembly histories of spirals is to determine whether the 
Local Group massive spirals are suitably typical by study- 
ing giant spirals beyond the Local Group. 

A number of surveys of the low surface brightness 
outskirts of spirals beyond the Local Group have been 
conducted recently. These surveys w ere however ei- 
ther sampling limited galactic volumes dMouhcine et al.l 
200 511 iBland-Hawthorn et alj|2005t iMouhcine et alj|2007t 
Ide Jong et al.ll2007D. or too shallo w to detect the old stel- 
lar tracers ( Davidgel 120061 120071 ). thus severely hamper- 
ing their impact. Measurements of galax y outskirts have 
been also attempted using integrated light (jMorrison et al.l 
1 19941 ) , and have succeeded on detecting low surface bright- 
ness tidal streams in the outskirts of a few n earby disk 
galaxies (jMartinez-Delgado et all 120081 I2009D . Those 
meas urements are however affected by large unce rtain- 
ties (|de Jong I [20081 ; iMartinez-Delgado et all 120081) . and 
are able to detect stellar structures down to a surface 
brightness limit of \i\ ~ 28 magarcseec -2 at the faintest 
(|Zheng et al.lll999l) . many magnitudes brighter than the 
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bulk of substr ucture detected in the outskirts of the Lo cal 
Group spirals (<Belokurov et alJl2007Hlbata et al.ll2007l ). 

Characterizing comprehensively the outskirts of spirals 
beyond the Local Group requires resolving panoramically 
the stellar content of those regions, giving access to the 
extremely low surface brightness structures. The required 
observations are however extremely challenging due to the 
photometric depth one has to reach, i.e., / ~ 26. — 28.5, 
to resolve old giant stars in galaxies beyond the Local 
Group, restricting this approach (using the present-day in- 
strumentation) to the small number of spirals closer than 
~10-12Mpc. 

Because the theoretical predictions are inherently sta- 
tistical in nature, and due to the stochastic nature of halo 
formation, constraining the assembly history of galaxy 
outskirts must rely in large part on measuring the de- 
mographics of their stellar populations. To this end, we 
have initiated a survey to resolve panoramically the stel- 
lar content of galaxy outskirts for a sample of nearby, i.e., 
up to ~ 10 Mpc, highly inclined spirals, distributed over 
a broad range of masses, i.e., circular velocities ranging 
from ~ 80kms _1 up to ~ 230 km s -1 , and morphologies, 
i.e., ranging from Sa to Sd. 

As part of this effort, we have targeted the edge-on 
galaxy NGC 891, often considered as the nearest prime 
analog of the Galaxy. Our group recently used deep op- 
tical imaging data, obtained with the Advanced Cam- 
era for Surveys on board the Hubble Space Telescope 
(HST/ACS), to investigate the properties of the resolved 
extra-planar stellar populations over the south-east quad- 
rant of NGC 891, extending up to ~ lOkpc from the 
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2009J). Succinctly sum- 
that the thick disk of 



NGC 891 shares comparable structural and chemical prop- 
erties to its Galactic counterpart. The stellar populations 
beyond the thick disk appear to possess significant small- 
scale variations in the stellar metallicity. Interestingly, 
those regions are found to be dominated by stars signif- 
icantly more chemically enriched than those populating 
the regions of comparable heights from the plane of the 
Galaxy. 

In the present contribution, we report the first results 
of our survey. Detailed analysis of the properties of the 
stellar content of different galactic components, the search 
for substructures, the determination of metallicity distri- 
bution functions and their spatial variation, and the prop- 
erties of the globular cluster system will be reported in 
forthcoming papers. The layout of this paper is as follow; 
in § [5] we present briefly the data set. In § [3] we report 
the discovery of a new morphological structures around 
NGC 891, and discuss the implications of our finding for 
the formation and evolution of spiral galaxies. 

Throughout the paper we use an intrinsic distance mod- 
ulus for NGC 891 of (m - M) Q = 29.94, with the tip of 
the red giant branch l ocated at I = 25.84 ± 0.04 ma g 
l|Mouhcine et al.ll2007t iTikhonov fc Galazutdinovall2005l) . 
NGC 891 is located at relatively low Galactic latitude 
(£ = 140.38°, b = -17.42°), and therefore it suffers from 
significant (though not la rge) extinction from foreground 
dust: E(B -V)= 0.065 (|Schlegel et al.lll998ft . 



2. DATA 

The Subaru Prime Focus Camera (Suprime-Cam) on 
the 8.2-m Subaru Telescope is a mosaic camera of ten- 
chip 2048 x 4096 charge-coupled devices, which covers a 
34 arcmin x 27 arcmin field of vie w with a pixel scale of 
0.20arcsec (jMivazaki et al.l 120021 ). The instrument was 
used to image the outer regions of NGC 891. In the fol- 
lowing, we will describe briefly the observations and the 
data reduction, the full details will be reported elsewhere. 

The observations were taken in the Johnson visual V- 
band and Gunn z-band. We have obtained a total of 10 
hours of good quality data in the V-band and 11.27 hours 
in the z-band, with seeing better than 0.6 arcsec, allowing 
us to resolve approximately the brightest two magnitudes 
of the red giant branch (RGB hereafter) of metal-poor stel- 
lar populations at the distance of NGC 891. By covering a 
vast ~ 90 kpc x 90 kpc region around NGC 891, the survey 
allows us to distinguish local density enhancements from 
large-scale structure of the halo of the galaxy and/or the 
foreground distribution of stars. 

The images were pre-processed using the the Cam- 
bridge Astronomical Survey Unit photometric pipeline. 
The final stacks were created by summing all pixels, 
weighted by the estimated seeing on each frame. The 
DAOPHOT / ALLSTAR software suite (|Stetsonlll987h was 
used to detect all sources down to 3cr above the sky. Iso- 
lated bright stars were identified over the frame to serve 
as point spread function (PSF) templates. The PSF was 
modelled as a Moffat function; experiments showed that 
allowing the PSF to vary spatially over the field did not 
improve significantly the fit to bright stars, so we adopted 
a spatially constant PSF. Finally, the instrumental magni- 
tudes were s hifted to agree with th e calibrated photometry 
presented in lReikuba et alj (|2009h . 

The analysis of the photometric errors indicates that 
the average i-band error varies from <jj « 0.05 at i » 
25.8 - 26.0, i.e., the bright end of the RGB, to er; w 0.15 
at i « 26.8 — 27.0, i.e., approximately a magnitude below 
the tip of the RGB, while the average color error varies 
from a(y—i) ~ 0.12 to <J(v~i) ~ 0.18 within the same 
magnitude range. To estimate the photometric complete- 
ness of the survey, we used the artificial star simulations 
performed for our analysis of the HST/A CS deep images 
of NG C 891, and described in full detail in lReikuba et al.l 
(2009). By comparing directly the number of stars de- 
tected in both HST/ACS and Supreme-Cam images as a 
function of magnitude and spatial density of RGB stars, it 
is possible to compute the photometric completeness of our 
survey. For stars located above ~ 6 kpc from the galactic 
plane, i.e ., regions where the c rowding is not at all impor- 
tant (see lReikuba et ail (|2009l ) for more details), the data 
are well above ~ 80 per cent complete for the upper one 
magnitude of the giant branch, i.e., i < 27. 

Objects were classified as artefacts, galaxies or stars 
according to their morphological structure on the im- 
ages. The point source catalog consists of NGC 891 RGB 
stars, NGC 891 asymptotic giant branch (AGB hereafter) 
stars, Galactic foreground stars, and unresolved back- 
ground galaxies. In order to isolate stars to probe the 
spatial distribution of the stellar populations in the out- 
skirts of NGC 891, we applied a series of magnitude and 
color cuts. RGB stars are defined as those with fore- 
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ground extinction-corrected i-band magnitude fainter than 
25.8 and brighter than 27, and with (V-i) color redder 
than 1.1, while AGB stars are selected as those with fore- 
ground extinction-corrected i-band magnitude ranging be- 
tween 25.0 and 25.8 and foreground extinction-corrected 
(V-i)o colors redder than 1.2. 

3. THE PANORAMIC LANDSCAPE AROUND NGC 891 

Fig. 1 shows the surface density map of RGB stars across 
the surveyed area. Numerous enhancements of the surface 
density of RGB stars are visible with the most striking 
being the large-scale complex of arcing loops and streams 
wrapped around the galaxy. About a half-dozen arc-like 
features are visible in the density map of old stars, ex- 
tending up to ~ 40kpc west and ~ 30kpc east of the 
galactic stellar disk. On the left side, corresponding to 
north-eastern part of the galaxy, it is possible to trace the 
full extension of a giant loop turning around, falling toward 
the plane of the galaxy, and extending further to the south. 
The spatial density of AGB stars does not exhibit however 
any obvious enhancements associated with the large-scale 
complex of streams of old stars. Note that there is no evi- 
dence for a significant clumping in the distribution in of ex- 
tended sources across our field, Galaxy clusters and other 
large-scale structure are unlikely to dominat e the counts at 
the a ngular scales of the observed features (|Caouch et al.l 
11993ft . 

It is tempting to suggest that the large-scale network 
of old star streams is connected, originating from a sin- 
gle accretion event, however it is impossible to argue this 
firmly based on the stellar density map alone. Neverthe- 
less, the arc-like features exhibit similar shapes, angular 
lengths, radii of curvature, separation from the parent 
galaxy, and levels of stellar density enhancements. In addi- 
tion, the stellar loops appear to be distributed and to cross 
each other in a pattern strikingly similar to the classical 
rosette-shaped loops predicted for tid ally disrupting dwarf 
galaxies in TV-b ody simulations (e.g. Ilbata fc Lewisl fi998; 
lLaw et al"1l2008ft . These properties provide circumstantial 
evidence that the arclike features around the galaxy may 
all come from a single accretion event (detailed modeling 
of the stellar stream will be presented elsewhere). The 
discovery of the NGC 891 giant stream in the first deep, 
panoramic survey of the Milky Way's nearest analog, to- 
gether with the previous discoveries of tidal streams in the 
Milky Way, Andromeda, NGC 5907, and NGC 4013, sug- 
gests that halo substructure in the form of tidal streams 
may be a generic property of massive spirals, and that the 
formation of galaxies continues at a moderate rate up to 
the present day. 

Fig. 1 shows that the spatial density of old stars varies 
along the streams, however nowhere among these fea- 
tures is there a region dense enough that it could be 
identified as the remaining core of the disrupted progen- 
itor dwarf satellite. The stellar streams are most likely 
the fossil remnants of a totally cannibalized system, al- 
though the core of the disrupted satellite could possi- 
bly be hidden behind the disk of the galaxy. Deep ra- 
dio observations of NGC 891 have reported the presence 
of a large gaseous filament and a number of counter- 
rotating clouds in its halo, extending to comparable dis- 
tances from the disk as the stellar stream (jOosterloo et al.l 



2007(). The structure and the kinematics of these gaseous 
structures appear to favor a scenario in which they are 
the results of a fiyby interaction with the gas-rich sate l- 
lite UGC 1807 (lOosterloo et al.ll200fl iMapeUi et al.ll2008ft . 
Fig. 2 shows the color-magnitude diagram of stars along 
the streams, with fiducials of RGB sequences of three 
Galac tic globular clusters from iDa Costa fc Armandroff 1 
(1990) superimposed. The morphology of the the color- 
magnitude diagram indicates clearly that the stream pro- 
genitor contained a broad range of stellar populations, 
with a mean metallicity comparable to that of 47 Tuc, 
i.e., [Fe/H] ~ —0.7, and ages older than a few Gyr with 
hardly any young stellar components, in sharp contrast 
wit h the typical stellar p opulation of gas-rich dwarf galax- 
ies (Skillman et al.ll2007|) . This indicates that it is highly 
unlikely that the large-scale complex of loops and streams 
in the halo of NGC 891 could be associated with the gas- 
rich dwarf companion. 

The second striking feature visible in the surface den- 
sity map of RGB stars is the extended super-thick envelope 
surrounding the galaxy. The defining high surface bright- 
ness components of a spiral, i.e., the bulge and the disk, are 
embedded in a vast flattened cocoon-like structure extend- 
ing vertically up to ~ 15 kpc, and radially up to ~ 40 kpc, 
with a structure quite unlike any classical notions of thick 
disk or halo. This structure is strikingly flatter than the in- 
ner component of the G alactic halo that dominates out to 
~ lO kpc above the disk (jPreston et al.lll99"ltlCarollo et al.l 

2008) . An additional highly flattened co mponent of the 
Gala ctic halo has been recently detected ([Morrison et al.l 

2009) . This halo flat component is however suspected to 
be con fined within very close distances from the Galactic 
plane (|Preston et al.lll99H iMorrison et~al~ll2009l ) 

A stellar structure resembling the Galaxy thick disk has 
been resolved in NGC 891, with a vertical scale height of 
hz ~ 1-45 kpc (jlbata et al.ll2009ft . Stars in this morpho- 
logical structure show no metallicity gradients both verti- 
cally and radially, and cover a wide range in metallicity, 
simil a r to what is observ ed for the Galaxy (jGilmore et al.l 
Il995t Hvezic et all 12008ft . However, at vertical distances 
larger than ~ 5 kpc from the plane o f the galaxy, star s 
belong to a more extended structure (|Ibata et al.l l2009h ■ 
and possess a negat i ve me tallicity gradient, albeit a mild 
one (|Reikuba et al.l [2009). The envelope surrounding 
NGC 891 is therefore highly unlikely an extension of the 
thick disk. Thus it appears that the flattened super-thick 
envelope surrounding the galaxy is a previously unknown 
component of NGC 891. 

Our previous investigations of the properties of the stel- 
lar content of the south-east quadrant, sampling stars of 
the super-thick envelope, find strong eviden ce for large 
amou nts of chemically-distinct sub-structures (jlbata et al.l 
2009|. This suggests that those regions are populated by 
numerous accretion remnants that are spread over a large 
volume and are still far from being fully phase mixed. 
Since it is natural to expect that the stellar populations 
in this quadrant are representative of those populating the 
super-thick structure, this indicates that the stellar struc- 
ture surrounding the galaxy was assembled from the tidal 
disruption of many accreted satellites. 

The super-thick structure presented here is not just a 
peculiarity of NGC 891: a preliminary analysis by our 
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group of deep panoramic observations of the nearby edge- 
on galaxy NGC 2683 shows clearly the presence of an 
almost identical structure surrounding the high surface 
brightness components of that galaxy. Detailed kincmat- 
ical measurements have revealed that a number of stel- 
lar substructures present at large radii of Andromeda co- 
rotate with the inner stellar disk, and are distributed in 
a gigantic flattened "extended disk" structure, extending 
radially from the high surface brightnes s inner disk o ut to 
~ 40kpc from the galaxy center (|Ibata et al.ll2005l ). i.e., 
strikingly similar to the radial extent of the super-thick 
stellar structure surrounding NGC 891. Given that we 
do not observe Andromeda edge-on, it is possible that the 
"extended disk" in that galaxy has a similar vertical extent 
to the new morphological structure identified in NGC 891 
and NGC 2683. 

The debris resulting from the disruption of a dwarf 
galaxy on a prograde orbit that is coplanar with the 
host galaxy d isk is expected to rela x into an extended 
rotating disk (jPenarrubia et al.l 120061 ). The detection of 
the Monoceros stream in the Galaxy ([Yannv et al.l [2003; 



llbata et al.ll2003D. the presence of an extended-disk around 
Andromeda ( Ibata et all 20051). and the stella r stream 
around NGC 4013 (|Martinez-Delgado et al.ll2009D . all sus- 
pected to be the remnants of the disruption of satellites 
moving in low inclined orbits, suggest that the accretion 
of satellites on low-inclination and low eccentricity orbits 
may be a relatively comm on process in the formation of 
spirals (jPohlen et al.l l2004i). Numerical simulations in the 
framework of the hierarchical cold dark formation scenario 
predict that large spirals experience several mergers of 
dwarf s with masses r anging from 1% to 10% of the host 
mass (IGaoet al.ll2004l). w ith many on low-eccentricity or- 
bits (jGhigna et al.lll998l) . The present work suggests that 
a super-thick stellar envelope formed by numerous accre- 
tions may be a common feature of large spirals. 
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Fig. 1. — Surface density map of RGB stars over the surveyed area around NGC 891. The over-densities of old RGB stars detected in 
the present study reveal a large complex of arcing streams that loops around the galaxy, tracing the remnants of an ancient accretion. The 
second spectacular morphological feature is the dark cocoon-like structure enveloping the high surface brightness disk and bulge. 
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Fig. 2. — i vs. (V — i)o color-magnitude diagrams for a stream field (left panel) and an off-stream field (right panel). Superimposed on each 
panel are RGB tracks (shifted to the distance of NGC 89 1) of three Galactic globular clu sters of different metallicities: M 15 ([Fe/H]=-2.2), 
NGC 1851 ([Fe/H]=-1.2), and 47 Tuc ([Fe/H]=-0.7) from Da Cost a fc Armandroffl (fl990h . 



